Canonical suppression in microscopic transport models 
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We demonstrate the occurrence of canonical suppression associated with the conservation of an 
U(l)-charge in current transport models. For this study a pion gas is simulated within two different 
transport approaches by incorporating inelastic and volume-limited collisions nn <-> KK for the 
production of kaon pairs. Both descriptions can dynamically account for the suppression in the 
yields of rare strange particles in a limited box, being in full accordance with a canonical statistical 
description. 
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I. INTRODUCTION AND MOTIVATION 

Since long, the various properties of excited nuclear 
matter occurring in heavy collisions have been studied 
within statistical descriptions 0, 0, 0, 0, IE • Here the 
underlying principle is the stringent assumption of full 
thermal equilibrium in the stage of the reaction to which 
the description is applied. 

In the last 10 years thermal model ansatzes of hadronic 
resonance gases via a grand canonical description have 
become very popular again because detailed experimen- 
tal results for the yields of individual hadronic species 
became available at various bombarding energies. In- 
deed, it was found that the statistical description can be 
applied very successfully for ultrarelativistic energies, in- 
dicating that the individual hadronic particles seem to 
evolve to a point of nearly perfect chemical equilibrium 
^| (f° r a comprehensive review the reader is referred 

However, for low collision energies (e.g. in the SIS, 
AGS and FAIR range) especially the strange particles do 
become rare. Thus, as has been well known for a long 
time, a canonical instead of a grand canonical description 
of the strange particles has to be applied. The canonical 
approach results in a suppression in the yields for rare 
particles in comparison to the grand canonical descrip- 
tion when treating the conservation of the corresponding 
U(l)-charge exactly [13,113,111 

Indeed, thermal models 
including the canonical suppression have been applied 
rather successfully for the description of the few mea- 
sured hadronic yields in heavy ion collisions at SIS ener- 
gies 

A dynamical interpretation of the canonical suppres- 
sion has recently also been offered by the formulation 
and solution of kinetic master-equations [l6l Il7j . It is 
the purpose of this investigation to show that such an 
occurrence of canonical suppression for rare particles is 
already warranted by present day transport models. 

A priori this is not a trivial statement. One could 
argue that the transport models are based on solving a set 
of coupled Boltzmann equations which originate from a 
grand canonical treatment. However, this is not the case 
as the realisations do conserve energy exactly and also 
the individual quantum charges like baryon number and 
(net) strangeness are conserved during the propagation 



and scattering processes as well. 

In the following section we briefly review the canoni- 
cal suppression of rare particles associated with the con- 
servation of an U(l)-charge and we also summarize its 
dynamical formulation via a master equation as given 
in |161 Il7| . In Section III we describe our simulation 
setup. A pion gas is simulated by incorporating inelastic 
and volume-limited collisions 7T7t <-> KK for the produc- 
tion of kaon pairs. We employ two completely differ- 
ent numerical realisations for the treatment of collisions. 
The first model under investigation is UrQMD 0, ^| , 
where two-particle collisions are realized via a standard 
geometrical interpretation. For previous studies of the 
thermodynamic properties of hadronic matter within the 
UrQMD model, the reader is referred to pol l2ll l23| . 
In the second model the collisions are treated by transi- 
tion rates in small spatial subcells. The latter algorithm 
has most recently been successfully introduced in a co- 
variant parton cascade to describe inelastic multiparticle 
Bremsstrahlung processes of type gg <-> ggg |24| . In Sec- 
tion IV we present the results of our analysis. Section 
V then provides a summary and a conclusion. We will 
argue that a (canonical) chemical equilibrium of kaons 
can by far not be achieved at intermediate energies for 
relativistic heavy ion collisions, at least if known cross 
sections for the production of kaons are assumed. 



II. CANONICAL SUPPRESSION AND ITS 
DYNAMICAL DESCRIPTION VIA A MASTER 
EQUATION 

The statistical description of systems incorporating the 
exact conservation of quantum numbers has been estab- 
lished for many years and there exist several approaches 
of varying complexity and generality [Tol llll . Here 
we give a brief summary of the results adapted to our 
needs. The conserved U(l)-charge to be considered is 
strangeness whose net value is taken to be zero through- 
out our calculations and simulations. Explicitly we con- 
sider only inelastic reactions of the type irir «-> KK, 
where the kaons and anti-kaons bear strangeness +1 and 
— 1 respectively. 

For a large number of kaons it is sufficient to treat 
strangeness conservation on the average (Nk) — (JVW) = 



2 



(with Nk, Nj^ being the number of kaons and anti-kaons, 
respectively). Introducing a chemical potential /i s , or fu- 
gacity X s = exp (^j, to control the net strangeness con- 
tent, one is led to the following grand canonical partition 
function 



Z* C (V,T, A s ) = exp (Zl 



x s z K + X s 1 Zj^j 



(1) 



where Z\ denotes the relativistic one-particle partition 
function for non-interacting particles of type i (pions or 
kaons in this case) , V is the volume, T is the temperature 
and rrii is the particle's mass. Ki denotes a modifed 
Bcsscl function. 



(2) 



In order to consider a more general case, the one- 
body partition functions in JJ| would have to be re- 
placed by sums over one-body partition functions of all 
hadronic particles carrying the corresponding integer in 
strangeness of 0, ±1, ±2 and ±3. 

A small number of kaons (roughly spoken when the 
average number becomes of order one) demands the 
conservation of strangeness to be treated exactly, i.e. 
Nk — Njr — for each state contributing to the partition 
sum. This constraint leads to a reduction of the available 
phase-space for the production process and ultimately 
one obtains the canonical partition function describing 
the system via (Iq and 1\ denote Bessel functions) 



Z c (V,T) = e X p(Zl)l Q (x) 



with 



2Z K . 



(3) 



(4) 



From Q and © one can calculate the density nx of 
particles with strangeness +1, the kaons, via the relation: 
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and considering the canonical case one finds 
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using that in the grand canonical picture (Nk) — (N-^) 



(6) 



requires A s = yJZyZ K = 1. 

Comparing (JSJ) and © one defines the canonical sup- 
pression factor < r] < 1 
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which contains all relevant information on differences 
in the equilibrium particle density between the grand 
canonical and canonical description. 

An alternative way of understandin g th is suppression 
is to consider kinetic master equations |l6lll7| by looking 
at a single process ab <-» cc. PN c (t) denotes the proba- 
bility to find N c particles c at a time t. In our case we 
have c = K and thus N c = N e holds exactly, whereas 
particles a and b (the 'pions') are assumed to be un- 
corrected and abundant. Furthermore, the probabilities 
for a single gain or loss process per unit time and vol- 
ume are denoted by G/V and L/V respectively, where 
G = (cg«) and L = (ulv) are the momentum- averaged 
cross sections for the gain and loss processes. With that, 
a master equation can be formulated |lfi| 



dP, 



G 



dt 



v (N a )(N b )P Na . 
!l(N a )(N b )P Ne + ^N?P Nc 



y(N c + l) 3 P Nc+1 



(8) 



A general solution of this equation is possible ^3] > but 
for our purpose it is sufficient to look at the kinetic equa- 
tion for the time evolution of the average number (N c ). 
It can be obtained when multiplying (jSJ by N c and sum- 
ming over it: 



d(N c 
dt 



(Na)(N b ) 



(Nl 



(9) 



This equation can be easily treated in the two limiting 
cases (N c ) ^> 1 and (N c ) -C 1. For abundant production 
of cc pairs, i.e. when (N c ) 1, the relation (iVf) ~ 
(N c ) 2 holds. On the other hand, for very rare production, 
i.e. when (N c ) <C 1, one has the relation (N%) ~ (N c ). 
Assuming a thermal momentum distribution with 
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and looking at stationary solutions for large times (the 
solutions describing an equilibrated system) one ends up 
with (e = G{N a )(N b )/L) 
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for (N c ) ^> 1. The opposite case (N c ) <C 1 leads to 
(N c ) 6 Z\Z\ _ c 
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(11) 



(12) 



Identifying particle-type c with kaons, it is clear that 
(|TT)l equals the grand canonical result J5J. (|T2"|) is just 
the leading term in an expansion of the canonical result 
©. Hence, it is verified that an abundant production of 
particles leads to a grand canonical description, whereas 
rare particle production requires a canonical description. 

Let us have a closer look at the canonical suppression 
factor (JJJ using the asymptotic behaviours 



lim 



h{x) 



x^oo Iq(x) 



1 



lim 



h{x) 



x^O Iq(x) 



X 

2' 



(13) 



3 



Here, one sees that in the grand canonical limit the 
kaon density is independent of the reaction volume, 
whereas in the canonical regime with the number of kaons 
(Nk) <C 1, it scales linearly with the volume as x oc V. 
Fig- H illustrates this behaviour. 



tion JHJl- Here the occurrence of a pair of rare parti- 
cles is explicitly taken care of. Due to the nature of 
pair production the existence of a kaon requires the exis- 
tence of an anti-kaon. This leads to an effective increase 
in the annihilation process compared to the Boltzmann 
Stosszahlansatz. Inspecting the loss term in @, the ap- 
proximation leading to a Boltzmann process is given by 




, n, Boltzmann process , 
(7V C 2 ) -» (N c 



(14) 



g.c. 



150 



200 



250 



300 



V [fm ] 



FIG. 1: Volume dependence of the kaon density for a canoni- 
cal ensemble (cf. eq. ©) for T — 170 MeV. g.c. denotes the 
grand canonical calculation, c. labels the canonical result. 



III. REAL TRANSPORT AND SIMULATION 
SETUP 



The main objective of our study is to demonstrate that 
current (large scale) transport models are able to repro- 
duce the effect of canonical suppression. This is not a 
trivial investigation. One might argue that the transport 
algorithms are based on strictly solving the Boltzmann 
equation which in turn originates from a grand canoni- 
cal treatment. In the stationary limit of the (coupled) 
Boltzmann equations the equilibrium phase space distri- 
butions, where the various collision terms for each par- 
ticle population vanish, obey Maxwell-Boltzmann statis- 
tics (or Bose-Einstein or Fermi-Dirac statistics when the 
Pauli factors are included), so that the individual densi- 
ties are that of a grand canonical ensemble. The reason 
for this is that for describing an annihilation within a 
Boltzmann description the two (rare) particles are as- 
sumed to be completely uncorrelated. This is however 
not the case in the description via the master equa- 



which then leads to the grand canonical description. 
Staying within the description of the master equation, 
the possibility for annihilation is much higher. When 
there exists a kaon, there has to be an anti-kaon with 
which an annihilation process is possible. (AT?) S> (N c ) 2 
for (N c ) <C 1, hence the effective annihilation rate incor- 
porated correctly in the master equation is much larger 
than for a Boltzmann description, leading to the canoni- 
cal suppression. 

On the other hand the various numerical realisations 
of the underlying transport equations do conserve energy 
and the individual quantum charges (like baryon number 
and net strangeness) exactly - as within a microcanonical 
treatment. The argument of an 'enhanced' annihilation 
should therefore also apply within these realisations. In 
order to address this question, we will concentrate on the 
volume dependencies in the grand canonical and canoni- 
cal regime respectively, as depicted in Fig. ^ 

The simulation setup consists of a large box of 20 fm 
side length holding a relativistic gas of pions. The pion 
gas provides a heat bath for a much smaller reaction vol- 
ume of variable size centered within the large box. In- 
side this smaller and likewise box-shaped reaction volume 
processes 

7T7T «-> K K 

are allowed, covering all possible isospin states of pions 
and kaons. The kaons are reflected by the walls of the 
small reaction volume and are thus bound to it, whereas 
these walls are permeable for the pions. An illustration 
of this special spatial setting is given in Figure El After 
equilibration, the kaon density within the reaction vol- 
ume should be governed by eq. (JHJ) which holds as a ref- 
erence for the transport model results. The idea therefore 
is to simulate different sizes of the inner reaction volume 
and to extract the number of kaons (Nk) for each simu- 
lation run by averaging over many timesteps in order to 
minimise the statistical fluctuations. These timesteps are 
sufficiently separated to avoid correlations. Furthermore 
the data are taken after the equilibration of the system. 

We implemented this scenario using two different types 
of transport descriptions - the microscopic transport 
model UrQMD [lg, yj| and a realisation of a stochas- 
tic transport model borrowed from a recently developed 
parton cascade [24|. The former model makes use of a 
geometrical interpretation of cross sections in order to 
solve the transport equations, whereas the latter relies 
on the explicit calculation of transition probabilities. 
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FIG. 2: Illustration of the two boxes for the numerical sim- 
ulation: The pions (diamonds) move inside the larger and 
fixed box and thus provide a heat bath for the kaons and 
anti-kaons (open and filled circles), which can only move and 
interact inside the small box. The size of the small box is 
varied for investigating the dynamical occurrence of canoni- 
cal suppression in small volumes. 



closely related to the formulation of the master equation 
© and its solution discussed in section |H] We have to 
stress one important point: The stochastic method is, in 
principle, flexible to introduce test particles, that is to 
'subdivide' each particle into a number N of testparticles 
|24j . However, for the following investigation it is crucial 
that the produced kaons are not subdivided into further 
testparticles. We will discuss this further below. 

The initial conditions in both schemes are chosen 
such that the pion gas acquires a temperature of T — 
170 MeV. The appropriate number of pions and the to- 
tal energy of the system are calculated via the use of a 
grand canonical partition function Q for pions alone, as 
kaons are absent in the initial state. The total energy 
evaluates to 

(16) 

A heat bath volume of 8000 fm , as used in our simu- 
lations, then corresponds to a population of 1348 pions 
bearing a total energy of 747.5 GeV. Initially each pion is 
assigned the same fraction of the total energy, giving one 
half of the particles momenta in the positive x-direction, 
while the remaining particles start out bearing momenta 
in the negative x-direction. The spatial distribution is 
random. 



IV. RESULTS FOR THE DYNAMICAL 
SUPPRESSION 



For standard applications, the UrQMD model provides 
full space time dynamics for hadrons and strings. It is a 
non-equilibrium model based on the covariant propaga- 
tion of hadrons and strings. All cross sections are fitted 
to available data or calculated by the principle of de- 
tailed balance. For our studies the code is modified such 
that only reactions tttt <-> KK , together with elastic col- 
lisions among the pions and kaons, remain possible. The 
7T7T — * KK inelastic reactions are assigned a constant 
cross section of either lmb, 5mb or lOmb, whereas the 
backward reactions are calculated via the principle of de- 
tailed balance. 

The simulation of 2 <-> 2 processes within the stochas- 
tic method is based on the calculation of a collision prob- 
ability for each possible particle pair per unit volume A 3 a; 
and unit time At via |24| 

At 

Here, v re i denotes the relative velocity and 022 is the 
cross section for the considered 2 <-> 2 process. Similar 
to the UrQMD setup, the cross sections are set to be con- 
stant in one direction and calculated via detailed balance 
for the reverse reaction. The so obtained probability is 
then compared with a random number between and 1 
to decide whether the collision does take place or not. 
The implementation of the stochastic model is therefore 



A necessary verification of the simulations reliability 
is to check for kinetic equilibration. Fig. [3] demonstrates 
that the energy distribution of the pions does become 
exponential over six orders of magnitude with the desired 
temperature of T = 170 MeV in both schemes. 

Figures 01 and now display the actual results in terms 
of the kaon density as a function of the (small) reaction 
volume. The minor fluctuations in the results indicate 
the small statistical errors due to the finite number of 
averages done. It is clear that the canonical suppression 
is reproduced by both transport schemes. The kaon yield 
is suppressed for small reaction volumes with respect to 
the grand canonical limit. This verification states the 
main result of the present study. 

In general, certain deviations from the theoretical val- 
ues are to be expected due to the choice of initial con- 
ditions. Pion number and energy content are fixed such 
that, without any kaons present, the system is tuned to 
the reference temperature of T = 170 MeV and fugacity 
A = 1. The production of the heavier kaons then leads 
to small changes in these parameters, as the total en- 
ergy in the system is conserved as well as the total num- 
ber of particles, because only particle number conserv- 
ing collision processes are included in the present study. 
This slightly affects the equilibrium number of pions and 
thus also the kaons. For the case of the stochastic cas- 
cade we investigated the relative deviations of temper- 
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FIG. 3: Logarithmic energy spectra for the pion gases within 
the UrQMD- and cascade-calculations. Both spectra are iden- 
tical within statistical fluctuations. 



ature T new = T + AT and effective mesonic fugacity 
X e f t = 1 + AA more closely. It is found that both devia- 
tions indeed grow with the number of kaons produced and 
thus with the size of the reaction volume. The observed 
deviations are on the order of two percent at most. Tak- 
ing these deviations into account, the theoretical results 
can be adapted to the actual conditions present in the 
system at a given size of the reaction volume by rescaling 
with the effective fugacity factor. The actual theoretical 
reference curves in Figs 0] and are corrected that way. 

The stronger deviations for the microscopical UrQMD 
model cannot be solely explained by the changes in tem- 
perature and fugacity. In fact the crucial effect when us- 
ing the geometrical concept of incorporating binary colli- 
sions is a decrease in the collision rate, when the expected 
mean free path A m .f. p . = (rao^) -1 for the particular re- 
action gets in the order of the interaction length */ ct 2 2 /7r 
as pointed out in [24| • This is due to the difference in the 
collision times of the involved particles viewed from the 
computational frame. During that interval the particle 
having the larger collision time must not collide again 
to ensure causality. Thus, the collision rate is decreased 
compared to the one given by the collision integral. As 
different densities are involved, one can not expect for- 
ward and reverse reactions to be affected in the same 
way and the changes in reaction rates lead to a shifted 
stationary 'equilibrium' value of the kaon density. 
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FIG. 4: Kaon density versus reaction volume as extracted 
from simulations within the stochastic cascade (triangles). 
For comparison, the dashed line indicates the grand canonical 
behaviour. The solid line shows the canonical volume depen- 
dence of the kaon density, based on © for a temperature of 
T = 170 MeV and corrected by taking small deviations in 
temperature and fugacity into account. 



This behaviour can be clearly seen in Fig. where 
the kaon density is depicted for different cross sections of 
the forward reaction irir — > K K. The smaller the chosen 
cross section and thus the larger the mean free path, the 
more accurately the theoretical result is reproduced. On 
the contrary, for typical real mesonic cross section the 
numerical shortcoming is in the range of 10 to 15 per- 
cent. Too few kaons relative to the pions are simulated in 
comparison to the theoretical limit either in the canonical 
or grand canonical regime. Such a discrepancy in kaon 
number relative to the pion number with respect to ex- 
perimental results has been reported in several URQMD 
calculations [2^. It might well be that a better imple- 
mentation of the collision criteria will enhance the kaon 
relative to the pion yield. We leave this for a detailed 
future investigation. 

Another instructive way to analyse the results is to 
look at the suppression factor 77 (7J as a function of the 
average number of kaons {Nk) c present in the system. 
This relation is independent of the temperature which 
can be easily shown: Since the suppression factor is a 
function of Z]^, r\ — r\{ZW the one-body partion func- 
tion can in principle be expressed in terms of 77, 

Z X k = Zjdv), (17) 
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FIG. 5: Kaon density versus reaction volume as extracted 
from simulations using the UrQMD model. Cross sections 
of lmb (triangles), 5mb (circles) and fOmb (diamonds) have 
been assigned to the forward direction of the inelastic reac- 
tions tttv <-> K K. The solid and dashed lines again provide 
the theoretically expected behaviour as seen in Fig. [I] 



which then leads to 



(N K y = v (n k )^ = nzjdn). 



(18) 



Fig. states the summary of the investigation and de- 
picts the above functional relation together with the sim- 
ulation results. Here the temperature and fugacity devi- 
ations have been taken into account for the cascade re- 
sults. For the calculations employing the UrQMD model 
the suppression factor has been calculated by taking the 
ratio of the actual kaon density nic(V) divided by the 
density tik(V — ► oo) in the limit of very large volumes. 
Within this prescription the sensitivity on the cross sec- 
tion vanishes almost completely. Apart from the demon- 
stration of the excellent agreement between theory and 
simulation, the plot illustrates the two limiting cases con- 
sidered when solving the kinetic equation @ . Small par- 
ticle numbers (N) <C 1 are associated with the canon- 
ical regime, whereas for abundant particle production 
(N) ^> 1 the grand canonical description is fully valid. 



V. SUMMARY AND CONCLUSIONS 

We have demonstrated that two current (large scale) 
transport models are able to reproduce the effect of 
canonical suppression. As a particular example we have 
chosen the production of kaons and anti-kaons via pions. 



FIG. 6: Suppression factor r\ as a function of the average kaon 
number. 



The overall kaon yield is suppressed with respect to the 
grand canonical limit for small reaction volumes and thus 
small numbers of produced kaons. 

The transport algorithms are based on solving the 
Boltzmann equation. Still they do conserve energy and 
net strangeness exactly. The canonical suppression for 
kaon numbers (iVV) considerably smaller than one, as 
pointed out in 16], originates dynamically from an en- 
hancement of the annihilation process by 1/(Nk) as 
compared to standard, grand canonical formulation of 
the Boltzmann equation. The reason is that any kaon 
in the system requires the existence of a correspond- 
ing anti-kaon due to strangeness conservation. Thus the 
probability of finding a particle anti-particle pair turns 
into a highly correlated conditional probability com- 
pared to that obtained via the uncorrelated Boltzmann 
Stosszahlansatz. The so enhanced annihilation probabil- 
ity then leads to the canonical suppression in the kaon 
yields. As the calculations show, it is manifest that the 
occurrence of canonical suppression is dynamically re- 
produced by the two presented schemes for solving the 
kinetic transport equations. A small discrepancy, i.e. a 
slightly lower kaon yield on the order of 10 to 15 per- 
cent, is present when employing typical cross sections. 
This discrepancy can be traced back to the numerical re- 
alisation of the occurrence of collisions via the standard 
geometrical description. As already emphasized, detailed 
investigations within the UrQMD model are in order to 
clarify how important such effects are for the understand- 
ing of kaon yields in relativistic heavy ion collisions, as 
typically too few kaons are produced relative to the pions, 
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when compared to experimental results |25|. 

Addressing especially the dynamical generation of 
canonical suppression, one has to worry about even more 
significant violations when other numerical transport 
schemes are invoked: One popular scheme is the so called 
test-particle method, i.e. each real particle is subdivided 
into N test-particles in order to obtain smoother distri- 
butions or to suppress other numerical artefacts (see e.g. 
[24^. When applying this idea of particle subdivision to 
the stochastic method, it will rescale the volumes V of 
Figure 01 by a factor 1/N, so that the grand canonical 
limit is achieved for much smaller volumes. This is, of 
course, unphysical. Thus, whenever the canonical sup- 
pression is of relevance for the understanding of particu- 
lar yields, the test-particle method has to be avoided. 

Another oftenly applied method is the so called per- 
turbative method with 'virtual' particles (for a discus- 
sion see e.g. H^). Here it is assumed that the exotic 
particles being produced are so rare that their behaviour 
does not alter the overall dynamics. If the backreaction 
of kaon production, i.e. the annihilation of a kaon and 
its antiparticle, is not considered correctly, then detailed 
balance is violated and the system can never achieve full 
and correct chemical equilibrium. 

When simulating real heavy ion collisions to address 
the production of kaons for lower or moderate relativistic 
energies (for a very recent work see |27j), one has indeed 
to worry whether these backreactions are important in 
the sense that the rates to achieve chemical equilibrium 
are smaller or comparable to the overall lifetime of the 
fireball. Various applications of thermal models for such 
low energies have been put forward in several works over 
the last years for descri bing hadronic yields in heavy ion 
collisions at SIS 

CI III 13 HI As only very few yields 
are available, the soundness of such an analysis is not 
evident. In principle, such an analysis rests on the as- 
sumption that chemical equilibrium among the individual 
hadronic particles is temporarily achieved, which implies 
that the corresponding time scales are sufficiently short 
compared to the lifetime of the fireball. 

From transport theory residing on binary elastic and 
inelastic hadronic collisions it has been known for a long 
time that for intermediate to moderate energies the kaons 
are predominantly produced before the initial motion is 
substantially degraded and when the system is still far 
from any (quasi-)equilibrium stage (for a discussion and 



review see e.g. 29] ). Putting it differently, when the mo- 
menta of the nucleons are sufficiently degraded and the 
system has to some extent thermalized, the timescale for 
production of strange particles via the considered inelas- 
tic kinetic reactions becomes extremely large. For a par- 
ticular calculation in a static environment at somewhat 
larger energies it was found that the equilibration of the 
kaons is exceeding the lifetime of a potential fireball by 
at least two orders of magnitude |3(J]. In this analysis 
the comparison was made, however, to a grand canonical 
estimate. 

On the other hand, as we will argue now, a chemical 
equilibrium situation of kaons with canonical suppression 
included cannot be achieved even temporarily when as- 
suming known cross sections. Within a grand canonical 
environment the rate for chemical equilibration of a par- 
ticle C for a reaction of type A + B ^ C + D is given by 
r -1 w (<TcDVcD)n<D, where no is the particle density of 
species D. With uk — n-j^ if there is only one kaon, then 
there also has to exist one anti-kaon in the particular re- 
alisation of the system. Hence, for such a situation of a 
canonical realisation, where the reactions 7r + 7r <-> K + K 
will drive the kaons to chemical equilibrium, the subse- 
quent equilibration rate is obtained by substituting for 
no the effective anti-kaon density 1/V. The rate thus 
reads 

T chcm:K ~ ( a KK ^it -n V K~k) y ■ (^) 

In an UrQMD simulation the typical averaged 
cross section rather close to threshold is given by 

( a KK-n™ V Klc) ~ 2mb ' S0 tliat T chcm;K ~ 5V^/(fm 2 c). 

Any typical reaction volume at intermediate time scales 
in a relativistic heavy ion reaction is in the order of 100 
fm 3 . Hence, kaons can not be expected to come or stay 
close to chemical equilibrium. 
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